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An analitical approximation of (0^) for a massive scalar field in a zero temperature vacuum state 
in static spherically symmetric spacetimes is obtained. The calculations are based on the method 
for computing vacuum expectations values for scalar fields in general static spherically symmetric 
spacetimes derived by Anderson, Hiscock and Samuel [Q,^ . The analitical approximation is used to 
compute {(f)^) in Schwarzschild and wormhole spacetimes. 
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. In the absence of a fully satisfactory theory of quantum gravity the semiclassical theory of gravity takes an important 
' role in studies of effects of the back reaction of the quantized fields upon the spacetime geometry. A primary technical 
^"5 difficulty in semiclassical gravity is that expectation values have the strong functional dependence on the metric tensor 
5^1^ and are generally impossible to be calculated analitically. For this reason, much efforts have been concentrated 
on developing approximate methods. 

In 1982, Page P] developed a method for obtaining an approximate expression for (Tp,y) for a conformal massless 
" ' ' scalar field in static Einstein spaces. Later a slightly different approach was proposed by Brown and Ottewill 0| 
I (see also Ref. p). In 1987, Frolov and ZePnikov constructed the approximate expression for (T)i^) for conformal 
ly-j ' massless fields in general static spacetimes; their construction was based on pure geometrical arguments and common 
(•~^ properties of the stress-energy tensor. The effective action approach had been elaborated by a number of authors 
M"|l5|- Recently, Anderson and Anderson, Hiscock and Samuel |^ developed a method of deriving an approximation 
' for (0^) and (T^i/) in a static spherically symmetric spacetime using the WKB approximation for the modes of the 
scalar field; the method is suitable for the fields with an arbitrary coupling to the scalar curvature; also the fields 
can be either in a zero temperature vacuum state or a nonzero temperature thermal state. 
O ■ The aim of this paper is to obtain the approximate expression for {(fP) using the Anderson-Hiskock-Samuel approach. 
^ '' Note that a practical applying of the method faces with considerable technical difficulty. It consists in necessity of 
. computing the mode sums and then expanding them in powers of e (where e — i{t — t') characterizes the point 
splitting). In Refs. such calculations was made approximately in the large uj limit. In this paper we present exact 
^ ' results on computing the mode sums. 

Another difficulty has no computational nature. The problem is that a zeroth-order solution, which determines an 
iterative procedure of solving the mode equation, can be only defined up to terms of the second order. Of course, 
^ ■ in principle, a final result of the iterative procedure should not depend on a particular choice of the zeroth-order 
■ ■ ' solution. However, the form of an approximate solution of the second or higher order should depend on this choice. 
In this paper we solve the mode equation iteratively using the zeroth-order solution which is slightly different from 
that in Refs. used. 

The paper is organized as follows. In Sec. II, following the Anderson-Hiscock-Samuel approach, we develop the 
second-order WKB approximation for an unrenormalized expression for (0^), where (/> is a massive scalar field with 
arbitrary curvature coupling in a general static spherically symmetric spacetime. In Sec. Ill the resulting expression 
is renormalized using the method of covariant point splitting. In Sees. IV and V we apply the resulting analitical 
approximation for (0^) to investigate the vacuum polarization of a massive scalar field in a Schwarzschild spacetime 
and in a wormhole spacetime, respectively. In Appendix A we derive the asymptotical expansions for the mode sums 
and integrals in the unrenormalized expression for (/>. In Appendix B some technical details are discussed. 

The units h = c — G = 1 are used throughout the paper. 
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II. AN UNRENORMALIZED SECOND-ORDER WKB APPROXIMATION FOR (c/.^) 



In this section we derive an unrenormalized expression for (0^) for a massive scalar field in an arbitrary static 
spherically symmetric spacetime. Note that at this stage our consideration will mainly follow the papers by Anderson 
and Anderson, Hiscock, and Samuel in which one may find some additional details. 

As in Refs. ||l|,^ , a Euclidean space approach is used. The metric for a general static spherically symmetric spacetime 
when continued analitically into Euclidean space is 



ds^ = f{r)dT^ + h{r)dr^ + r^{de'^ 



}dip'^). 



(1) 



Here r = it is the Euclidean time, and / and h are arbitrary functions of r which, if the space is asymptotically flat, 
become constant in the limit r —^ oo. 

Consider a quantized scalar field 4> with mass m and coupling ^ to the scalar curvature R. We assume that the 
field is in a vacuum state defined with respect to the Killing vector which always exists in a static spacetime. An 
unrenormalized expression for {(fp) can be computed from the Euclidean Green's hmction Ge(x,x). This expression 
given in Refs. |l|,|^ is 



((/)^(a;))unroii = Ge{x, t; x, f) 
1 

= / duJCOs\Lj(T 



1=0 



(2/ + l)puiquji 



where the modes p^^i and q^ii obey the equation 
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They also satisfy the Wronskian condition 
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(3) 
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where C^ji is a normalization constant. 

Let us stress that points in Eq. (|^) are splitted. Namely, points are separated in time so that e = [t — f ), r — r, 
6 = 6, (p = ip. As was first pointed out by Candelas and Howard |T^] for the case of Schwarzschild spacetime, the 
Euclidean Green's function have superficial divergences with this separation of points. As discussed in Refs. 
these cannot be real divergences because the Green's function must be finite when the points are separated; to remove 
the divergences one has to subtract some additional counterterms. The terms (r/^/^)^^ in brackets in Eq. (j^) are 
such counterterms (for details, see discussion in Refs. 

There is a WKB representation for the modes which is very useful in calculations of {(jy^)- The WKB representation 
is obtained by the change of variables 



P-' = 77e^TFu7iexp 



(2r2VF)i/2 



1/2 



dr 



{2r^W) 
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(5) 



dr 



where W is a new function of r. Substitution of Eq. (||) into Eq. (^) shows that the Wronskian condition is obeyed 
if Cuji = 1, and substitution of Eq. (|^) into Eq. gives the following equation for W: 



f 



where the scalar curvature R is 
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and the prime denotes the derivative with respect to r. The usual way to work out Eq. (|^) is to solve it iteratively. 
It is obvious that a choice of zeroth-order iteration determines the further iterative procedure. In Refs. it was 
suggested to take 



uj^ +m^f + {1 + 1/2)' 



f 



1/2 



as the zeroth-order solution. This choice seems to be convenient because it simplifies some further calculations. 
However, this is not the best choice. Really, consider the asymptotically flat region of spacetime where f,h^l and 
/', /", h', h" 0. Assume also that W' , W" there. In this case the equation ra) reduces to 



with 



1/2 



Hence, it is seen that a natural choice for the zeroth-order solution is ly^") — Q. Note that W'^^^ — 
So, we take = 17 as the zeroth-order solution. Then, the second-order solution of Eq. (ptj is 



(7) 



w = n + = n - 
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with 
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Discuss the question about applicability of the second-order solution. It is correct provided ^ This 
condition must be fulfilled for each mode. In particular, for the lowest null mode with uj = and I = we obtain 
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(9) 



In the other words, the last condition means that the Compton length Tc = m^^ is much smaller than the characteristic 
length scale of the spacetime geometry. Restrictions, which the inequality (|^) imposes for the metric functions / and h, 
ensure correctness of the WKB approximation. Hereinafter we shall suppose that / and h answer these requirements. 

To obtain the second-order WKB approximation for {(/)^) we substitute Eq. (||) into (||), and then into (H). Neglecting 
terms of the fourth order and higher we can finally find 
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11 + 1/2 
8 r^n^ 
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with 



U{r) ^ m^f + l{l + l)- 



III. A RENORMALIZED EXPRESSION FOR {(j,'^) 

In this section we derive a renormalized expression for (0^) in the framework of the second-order WKB approxi- 
mation. 

As is usual for the method of point splitting, in order to renormalize ((^^)unron one should subtruct renormalization 
counterterms from ((/)^)unron and take the limit e (r' r). Schematically, 
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The renormalization countertcrms for a scalar field has been obtained by Christensen which used the DeWitt- 
Schwinger expansion for the Feynman Green's function. In the case of a massive scalar field they are given by 
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Here a is equal to one half the square of the distance between the points x and x' along the shortest geodesic connecting 
them. C is Euler's constant, R^p is the Ricci tensor and ct" = ct'". 

To perform the procedure of renormalization in practice, one should expand the expressions for (02)unrcn and (02)ds 
in powers of e. The expansion of (02)ds can be easily found by using the expansion of a and its derivatives in powers 



of e given in Refs. ||l|. 



a=-/6^ + 0(64), a 



24fh 



f/2 

e3 + 0(e5). 



(13) 



In general, we have (02)ds = a€~^ + 61ne + c + 0{e^) where first two terms describe ultraviolet divergences. It 
is well-known that the second-order WKB approximation for (02)unicn contains the same divergences, and so the 
subtraction of {4>'^)-ds annihilates terms of (</'2)unrcn that is infinite in the limit e — > 0. 

To obtain the expansion of {(jy^) um-cn in powers of e we rewrite Eq. ( |l0[ ) in the form being more convenient for 
further analysis: 



with 



47r 



2 / J,2\ 
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—So{e,fi) - 5*1(6, /i) — — S'2(e,Ai)j 



f f f 



1 ffV f fV 
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fV f f 



y/u^+fi^ + {l + l/2y 



1 ff\' f f 



(14) 



(15.1) 



'S'„(e,/i) = 



n = 1,2,3, 



du cos 



(/ + 1/2)2 



^ K+;i2 + (; + l/2)2]»+l/2' 



(15.2) 



N-ie,,)^l ducos(u'-^)t,-,,-^^ (15.3) 



n= 1,2,3,..., 
m = 1, 3, . . . , 2n — 3 



where we denote 



M2 = m2r2-i. (16) 



Note that ^? >{)\ir> (2m)" 
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Now one must calculate the asymptotical expansion for the functions So{e,ii), Snie,fj,) and N^{e,fi) in the limit 
e — > 0. In Appendix A details of such calculations are given. Here we present some final formulas: 



^) - - (2n -l)!! r y + ^) + (2;H]TI + 1^)' 



where (2n — 1)!! = 1 • 3 • 5 • • • (2n — 1), n = 1, 2, 3, . . . and So(2nfi) and 5„(27r^) are determined by Eqs. ( |A.14D . 
Note that the functions N^{e, fj.) are regular in the limit e ^ 0. Hence one may directly put e = 0, so that 

Changing the order of summation and integration in Eq. ( p^ and integrating over u gives, in particular, 

9 9 °° / + i 



15 



(19.2) 



where i = 1, 3. 

Note also that the functions 5„(27r/i) and N^{^) have a simple asymptotical form at large values of /x, or, corre- 
spondingly, at large values of mr (see Appendix B for details): 
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Substituting Eqs. ([l7|,|l9.l|) into Eq. (|l|) we may find the asymptotical expansion for ((/> ) 
e ^ 0. Finally, earring out the procedure of renormalization (^Tj), i.e. subtracting (0^)dSi wc obtain the renormalized 
expression for (0^) in the framework of the second-order WKB approximation: 
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IV. {(j)'^) IN A SCHWARZSCHILD SPACETIME 



In this section we apply the analitical approximation for {(jP'), obtained above, to investigate the vacuum polarization 
of a massive scalar field with arbitrary curvature coupling in a Schwarzschild spacetime. 
For a Schwarzschild spacetime the metric functions / and h are 

/ = k-^l-^. (22) 

r 

where M is the mass of the black hole. The event horizon is at 

Tg = 2M. (23) 

The scalar curvature R is identically zero for a Schwarzschild spacetime. 

We shall compute {(fP') in the region exterior to the event horizon where the spacetime is static, i.e. r > r^. The 
condition ensuring applicability of the WKB approximation in this region, can be rewritten as 



2 

Is' 



» 



1 1 

V ^ 32p4(i 



(24) 



where Tc = is Compton radius corresponding to a scalar field with the mass m, and p = r/vg is the dimensionless 
coordinate, p > 1. For example, let us take rg = 2Mq = 1.48 • lO^cm and Tc = = 3.86 • 10~^^cm where Mq is 
the Sun mass and nic is the electron mass. In this case the inequality (24) is fulfilled for p — 1 > 10^'^". 

We shall assume that rg/vc » 1, and so mr = {rg/rc)p » 1 because of p > 1. Therefore, the functions Sn and 
iV™ can be taken in their asymptotical form (EQ) and the expression (p|) for {<jp) can be written down as follows: 



4 V 4A2p2 J I6p2 1920 A2p4 48 A2p5 7680 A^pS 

«i Z^-, 13 , 35 5 \ , a2 A 3 9 , 

1 - - + TT ' (25) 



AV(1-P~^)^ V 3p 6p2 2pV A6p^(l-P"^) V P 4p 

where A = 2Mm = Vg/rc, ai = 3, 646 • 10^"^ and 02 — 1.281 • lO^'^. Note also that at large values of p, i.e. far from the 
event horizon, the last expression has a simple asymptotical form. Taking into account that ln(l — i) — (i + 
for large values of x we obtain 

16.^M2(0^)„n«-^^^. (26) 
Some results of numerical calculations using the analitical approximation (|2^) are shown in Fig. 1. 



V. (02) IN A WORMHOLE SPACETIME 

In this section we apply the analitical approximation (|2l] ) to calculate (02) for a massive scalar field in a wormhole 
spacetime. 

Here we shall regard a wormhole as a time independent, nonrotating, and spherically symmetric bridge connecting 
two asymptotically flat regions. The metric of wormhole spacetime (continued analitically into Euclidean space) can 
be taken in the form, suggested by Morris and Thorne |l8[ : 

ds^ = f{l)dT^ + dl^ + r^{l){d0'^ + siiP 9 difi"^), (27) 

where I is the proper radial distance, / e (— cx),+oo). We assume that the redshift function f{l) is everywhere finite 
(no event horizons); the shape function r{l) has the global minimum at Z = 0, so that vq = imn{r{l)} = r(0) is the 
radius of the wormhole throat. In order for the spacetime geometry to tend to an appropriate asymptotically flat 
limit at Z ^ ±00 we impose 

lim {r{l)/\l\} = 1, and lim f{l) = /±. (28) 

/ — >±00 I — *-±CXD 

For simplicity we also assume symmetry under interchange of the two asymptotically flat regions, / ^ — /, that is, 
r(0=r(-/) and/(0 = /(-0- 
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Note that two metrics, (|l|) and (27), can be derived from each other by substitution 



Now 



dr _^ 1 



1 dh 



h dr 



(29) 



(30) 



where the prime denotes the derivative with respect to I. 

Let (f) he a, massive scalar field with minimal curvature coupling (^ = 0) in the wormhole spacetime. Consider the 
simplest example, given in Rcf. p8|, when the metric functions f{l) and r{l) are 



/(o^i, r(o=y^ 



(31) 



The condition (^, ensuring applicability of the WKB approximation in the wormhole spacetime, reduces in this case 
to 



1 r 



1 



ra - — = - 



2r 2(/2 + r2)- 



(32) 



It is seen that this inequality is satisfied for all values of / provided r\ ^ {2rn)~^ = rc/2. This means that the throat's 
radius has to be much greater than Compton radius. If so, then mr ^ 1, and the functions Sn and can be taken 
in their asymptotical form (Eff). Now the analitical approximation (pl^) for (0^) can be written as 



167rV, 



2^2/ i2\ 
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Pi 



,2^3 
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63 H 



mii-l 1 



1 1 



Ap-^+pl) ' ^P^+pI (p2+p2)2 (p2+p2)3 



P^+pI ' (p2+p2)2 



(33) 



where p = ml the dimensionless proper radial distance, po — mro is the dimensionless throat's radius, and the 
numbers bi are h = 62 = 8.51 • 10~3, 63 = i|, 64 = 5.851 • 10"^, 65 = 5.126 • 10"^. Note also that far from the 
wormhole's throat, p — > 00, the analitical approximation ( |3^ ) has the simple asymptotical form: 
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2^2 li2\ 
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(34) 



which does not depend on the throat's radius. 

In Fig. 2 we present some results of numerical calculations obtained by using the analitical approximation (|33|). 
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APPENDIX A: ASYMPTOTICAL EXPANSION FOR 5'jv(e,/i) IN THE LIMIT e ^ 



In this appendix we derive the asymptotical expansion for the functions So{e, p) and S',i(e, p) in the limit e ^ 0. 
Consider the following sums over I: 



So 



1 + 1/2 



+ + 1/2)2 



(A.l) 



(1 + 1/2) 



2n-l 



^ [y2 + (/ + l/2)2]n+l/2' 



(A.2) 
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where n = 1,2,3, .. . and y > 0. To transform them we use the Abel-Plana method of summarizing of convergent 
series The essence of this method consists of following: the series is presented as a contour integral 



res 



/(z)cot7r(z- -) 



/(z)cot7r(z - -)dz, 



(A.3) 



where C is a contour which surrounds a region on the complex plane containing only the poles of cot7r(z — the 
poles of f{z) (if they exist) have to lie outside of the contour C. Deforming the contour by a special way and currying 
out not difficult algebraic transformations gives the following relation: 



^/a+i)=limj/ f{x)dx + 



1^0 



f{-tt)dt 



f{tt)dt 



Using this formula we may transform Eqs. (AJ) and (A^) as follows: 



so(y) 



dx + 2 



dt 



v/y^^^ +1 



(A.4) 
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where 



(2/2 +2;2)"+l/ 



;iTT7i-2(-i)" 



j.2n-l 



dt 



(y2 _ t2)n+l/2 e27rt + 1 ' 



Calculating the first integral in Eq. ( [A.5| ) and earring out an integration by parts in Eq. (A. 6) yields 

ry t dt 



so{y) = -y + 2 



(A.6) 



(A.7) 



sn{y) 



y dt d fl dV'-^ f i2("-l) 



(2n-l)!!yo y/y^^dt \ t dt 



The expressions for the functions S'„(e,/i), given by Eqs. (15.1) and (15.2), may be rewritten as 



Snie,IJ,)^ / duCOs{us)Sn{\/ u2 + ^2)^ 



(A., 



where e = er~^f^^'^ and n = 0, 1, 2, . . .. 

In order to derive the asymptotical expansion for So{e,iJ,) and >5'„(e,/i) in the limit e — > we use the following 
asymptotical relations: 



/i(e) = / ^u2 + ^2 cos{ue)du 
Jo 



1 ^' efi 



+ Oie^ Ine), 



(A.9) 



/2(e) = / du cos(ue) / 
Jo Jo 



f{t)dt 



poo 1 pOO 

To prove the first formula (|A.9|) we make the following transformations: 



(A.IO) 
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(w + /i ) cos(ue) 
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y/yF+J^ 

2 \ I ^ cos(ue) 
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du 



Here Kn{y) are the Bessel functions of the second kind. Taking into account the asymptotic of Ki{y) in the hmit 
2/^0: 



y 2\ 2 2 



O(y^lny), 



we obtain the relation 

To prove the second formula (A. 10) we change the order of integration over u and t in I2 and compute the integral 
over u, so that 



cos(u£) du 



cos(ue) du 



pfj. poo 

h{e)= f{t)dt 7:"^^'r + / !{^dt — 



(A.11) 



Here Ko{y) and 1^0(2/) arc the Bcssel functions of the second kind. Taking into account the asymptotics of Kq and Yo 
in the limit y 0: 

^0(2/) = -(ln| + C) + 0(2/2 In y), 
Yo{y)^-{\n^ + C) + Oiy'lny), 



we can obtain the relation (A. 10) 



Finally, substituting Eqs. (|A.7[ ) and (|A.8|) into Eq. ( |A.9D and using formulas (|A.9|) and ( |A.10D gives 



5o(e, = 7; + ^ - ^) f + Cj - Y - A*'^o(2^A*) + 0(6^ In e), 



(A.12) 
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(2n-l)!! V 2 / (2n-l)!! 
where {2n — 1)!! = 1 • 3 • 5 • • • {2n — 1), 71 = 1, 2, 3, . . . and 5o(27r/i) and 5„(27r/i) are 
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APPENDIX B: ASYMPTOTICS FOR S'iv(27r^i) AND Afji^^(M) 



In this appendix we obtain asymptoti cs at l arge values of argument for the functions S'„(27m]and N^(fi). Consider 
the functions S'„(27r/x) defined by Eqs. (A.14). Making the substitution y = 2itiix in Eqs. (A.14) we find 
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yln|l - A2y2l 
ev + 1 



1 d y-' 



(B.l) 



ev + 1 



where A = (27ry^)~^ and n = (rn^r^ + ■ Note that the integrands in Eqs. (B.l) contain exponential functions 
and are exponentially decreasing at large values of y. Hence the main contribution into the integrals is provided with 
values of integrands in the region < y < 1. We are interesting in the case ^ 1. In this case A « 1 and Ay <K 1 if 
< y < 1. No w we may use the asymptotical formula ln(l — A^y^) = — A^y^ — ^A^y"* — |A^y^ + 0(A*y®). Substituting 
this into Eqs. (B.l) and integrating gives 

7 ^ . 1 
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52(2^^^) = ^- + 0((mr)-6), S^{2ti^i) = + 0((mr)-« 



+ 0((mr)-4), 



(rnr)'* 

Consider now the functions Nl'^{fj,): 

oo 
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where i = 1,3. Using the formula (5.1.26.23) from Ref. ||2^ we find 
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(B.3) 



V/ ( i - z/i ) - ij' 
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where = r'(z)/r(z) is the digamma function. It is also easily to see that 

and Nli^,)=NU^,)-f^^N^i^,). 

Using the asymptotical properties of the digamma function (see, for example, Ref. [^) we can obtain for large values 
of fi the following asymptotics: 



1 



2(mr)2 



+ 0((mr)-4), 
+ 0((mr)-*). 



1 



4(mr) 



0((mr) 



(B.5) 



P. R. Anderson, Phys. Rev. D 41, 1152 (1990). 

P. R. Anderson, W.A. Hiscock, and D.A. Samuel, Phys. Rev. D 51, 4337 (1995). 
D. N. Page, Phys. Rev. D 25, 1499 (1982). 

M. R. Brown and A. C. Ottewill, Phys. Rev. D 31, 2514 (1985). 

M. R. Brown, A. C. Ottewill, and D. N. Page, Phys. Rev. D 33, 2840 (1986). 

V. P. Frolov and A. I. Zel'mkov, Phys. Rev. D 35, 3031 (1987). 

V. P. Frolov and A. I. Zel'nikov, Phys. Lett. B 115, 372 (1982). 

V. P. Frolov and A. I. Zel'nikov, Phys. Lett. B 123, 197 (1983). 

V. P. Frolov and A. I. Zel'nikov, Phys. Rev. D 29, 1057 (1984). 

V. P. Frolov, Proceedings of the Lebed ev Institute of the Academy of Science of USSR 169 (1986). 



V. Frolov, P. Su tton and A. Zel 'nikov, ^r-qc/9909086 
I. G. Avramidi, ^iep-th/951o"l4c| 



I. G. Avramidi, Teor. Mat. Fiz., 79, 219 (1989). 



10 



[14] I. G. Avramid i. Nucl. Phys B 355, 712 (1991). 
[15] J. Matyjasek, ^r-qc/9912o"2C . 

[16] P. Candelas and K.W. Howard, Phys. Rev. D 29, 1618 (1984). 

[17] S.M. Christensen, Phys. Rev. D 14, 2490 (1976). 

[18] M.S. Morris and K.S. Thorne, Am. J. Phys. 56, 395 (1988). 

[19] M.A. Evgrafov, Analytic Functions (Nauka, Moscow, 1968) [in Russian]. 

[20] A. P. Prudnikov, Yu.A. Brychkov, and O.I. Marichev, Integrals and Series, voL 1, (Gordon and Breach, Amsterdam, 1996). 
[21] M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions, (US National Bureau of Standards, Washington, 
1964). 



11 



-5e-05 - 



-0.0001 - 



-0.00015 - 




FIG. 1. The curves in this figure display the value of (0^) for a massive scalar field in the Schwarzschild spacetime. From 
bottom to top the curves correspond to the values A = r^/rc = 10, 20, 30. 
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FIG. 2. The curves in this figure display the value of 
to top the curves correspond to the values po = ro/rc = 



0^) for a massive scalar field in the wormhole spacetime. From bottom 
10, 12, 14. 
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